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ABSTRACT 

The Type lib supernova (SN) 1993J is one of only a few stripped-envelope SNe with a progenitor star identified in 
pre-explosion images. SN lib models typically invoke H envelope stripping by mass transfer in a binary system. 
For the case of SN 1993J, the models suggest that the companion grew to 22 M 0 and became a source of ultraviolet 
(UV) excess. Located in M81, at a distance of only 3.6 Mpc, SN 1993J offers one of the best opportunities to 
detect the putative companion and test the progenitor model. Previously published near-UV spectra in 2004 showed 
evidence for absorption lines consistent with a hot (B2 la) star, but the field was crowded and dominated by flux 
from the SN. Here we present Hubble Space Telescope Cosmic Origins Spectrograph and Wide-Field Camera 3 
observations of SN 1993J from 2012, at which point the flux from the SN had faded sufficiently to potentially 
measure the UV continuum properties from the putative companion. The resulting UV spectrum is consistent with 
contributions from both a hot B star and the SN, although we cannot rule out line-of-sight coincidences. 
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1. INTRODUCTION 

The Type lib supernova (SN lib; see Filippenko 1997 for a 
review) subclass exhibits spectroscopic evidence for hydrogen 
emission at early times (the defining property of Type II 
supernovae, SNe II), Hei absorption at later times (as in 
SNe lb), and renewed broad Ha emission in the nebular phase 
(e.g., Filippenko 1988; Matheson et al. 2000; Taubenberger 
et al. 2011; Shivvers et al. 2013). While SNe lib require 
most of the progenitor’s original H envelope to be lost prior 
to the SN explosion, a low-mass (0.01-0.5 M Q ) H envelope 
must be retained. Single massive stars (e.g., >30 M e ) could 
potentially lose their envelopes in a stellar wind, but models 
indicate that only a limited range of initial masses would result 
in a very low-mass envelope at the time of the explosion 
(Podsiadlowski et al. 1993; Woosley et al. 1994 and references 
therein). More likely progenitor scenarios involve lower-mass 
stars (e.g., evolved red supergiants) that lose their envelopes 
during mass transfer to a binary companion (Podsiadlowski 
et al. 1993; Nomoto et al. 1993; Woosley et al. 1994). The 
binary progenitor scenario is also consistent with SN rates, 
which strongly suggest that a substantial fraction of stripped- 
envelope SNe, including SNe lib, must originate with relatively 
lower-mass stars (<25 M 0 initially) stripped in binary systems 
(Smith et al. 201 1). 

Along with SNe 2013df (Van Dyk et al. 2014), 201 ldh 
(Maund et al. 2011; Van Dyk et al. 2011, 2013; Ergon et al. 
2014), and 2008ax (Li 2008; Crockett et al. 2008), SN 1993J 
is one of just four SNe lib with progenitors directly identified 


10 Deceased 201 1 December 12. 


in pre-explosion images (Aldering et al. 1994; Cohen et al. 
1995). The progenitor star, a K supergiant, was consistent with 
SN lib models that invoke supergiant mass loss onto a binary 
companion. By 2006, the SN had faded sufficiently for Maund & 
Smartt (2009) to confirm that the K-supergiant signature in the 
spectral energy distribution (SED) had disappeared. Even before 
this disappearance, however, excess flux in the near-ultraviolet 
(NUV) and B bands suggested the possible presence of the hot 
companion, thereby supporting the binary-star model (e.g., Van 
Dyk et al. 2002; Maund et al. 2004). 

However, the specific source of the NUV excess was not im- 
mediately obvious. SN 1993J has a number of stars within 1" of 
it (Figure 1). Using Hubble Space Telescope ( HST ) images, Van 
Dyk et al. (2002) initially showed that some, but not all, of the 
excess NUV could be explained by previously unresolved stars. 
Maund et al. (2004) later obtained an optical/NUV spectrum of 
SN 1993J and the putative companion with the Low Resolution 
Imaging Spectrometer (LRIS; Oke et al. 1995) on the Keck-1 
telescope and detected spectral absorption features consistent 
with a massive B-type supergiant. Even with a (y.'7-wide slit, 
however, these spectra would have been contaminated by Stars 
E, F, G, H, and I (Figure 1). Maund et al. (2004) calculated an 
upper limit on the contribution from the surrounding stars by 
simulating the spectral profile and concluded that the spectral 
features are dominated by a source at the position of the SN. 
Line-of-sight coincidences were not ruled out. 

While these results further support the binary scenario for 
SN 1993J, a direct detection of continuum from a hot B-star 
companion remains unconfirmed to date. Maund & Smartt 
(2009) predicted that with the observed rate of SED decline, a B 
star at the location of SN 1993J could be detected in the U and B 


1 


The Astrophysical Journal, 790:17 (13pp), 2014 July 20 


Fox ET AL. 


SN 1993J 
(a) WFC3/F336W 




-K 



J 


N 

A 


■ A 


pE 



SB 




0.5 

l 


H 


L 


t 


SN 1993J 
(b) WFC3/F555W 



! SN 1993J 
(c) WFC3/F814W 



SN 1993J 
(d) COS/MIRRORA 







Figure 1. Field of SN 1993J imaged in (a) WFC3/F336W, (b) WFC3/F555W, (c) WFC3/F814W, and (d) COS NUV, all to approximately the same scale and 
orientation. The dashed black circles correspond to the COS 2" 5 diameter aperture, which can also be considered poor seeing from the ground. Star labels designate 
the SN and other stars previously identified in HST images by Van Dyk et al. (2002) and Maund et al. (2004), together with additional fainter stars. 


bands by 2012. Here we present / 1ST U V, optical, and infrared 
(IR) photometry and spectroscopy of SN 1993J obtained in 
201 1 November through 2012 April. Section 2 summarizes the 
optical/IR photometry, Section 3 presents our Keck spectrum, 
and Section 4 describes the UV spectroscopy and spectral 
extraction procedures. We analyze the data in Section 5 with 
B-star models. Our results and conclusions are discussed in 
Section 6. 

2. WFC3 OBSERVATIONS 

SN 1993J was observed with the HST Wide-Field Camera 
3 (WFC3) UVIS and IR channels in 2011 November-2012 
February as part of program GO-12531 (PI: A. Filippenko), 
as summarized in Table 1. Figures l(a)-(c) display the WFC3/ 
F336W, F555W, and F814W images. Star labels designate the 
SN and other stars previously identified in HSTV images by 
Van Dyk et al. (2002) and Maund et al. (2004), together with 
additional fainter stars. 

Photometry was extracted from the individual WFC3 “fit” 
images in all bands using Dolphot v2.0 (Dolphin 2000). The 


input parameters are those recommended by the Dolphot WFC3 
Users’ Manual. Aperture corrections were applied. The F438W, 
F555W, and F625W data, all obtained on 201 1 December 24 in 
a set of orbits and therefore with the same orientation, were 
processed together. All other bands were obtained at other 
orientations and hence were processed in separate Dolphot runs. 
The data acquired on 2012 February 19 in the F814W and 
F850LP bands exhibit residuals created by a bright (red) source 
or sources likely observed immediately prior. The backgrounds 
around the SN environment in each of these bands are slightly 
elevated relative to the overall backgrounds in these images. 
Since Dolphot was set to measure the sky background locally, 
we do not expect these elevated backgrounds to significantly 
affect the resulting photometry. The F336W data obtained in 
the same sequence on that date appear not to show residuals. 
The resulting magnitudes in the WFC3 flight system (Vegamag) 
are listed in Tables 2 and 3 and plotted in Figure 2. 

We also detect the presence of expanding light echoes in the 
field of SN 1993J. Figure 3 shows F555W filtered images of 
SN 1993J taken in both 2002 and 2012. The 2002 data were 
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Table 1 

HST GO- 12531 COS/WFC3 Photometric Observations 


UT Date 

Epoch 

(days) 

Instrument 

Filter 

Central Wavelength 

(A) 

Exposure 

(s) 

20120331 

6944 

cos 

MIRRORA/G140L 

1105 

7466 

20120406 

6950 

cos 

MIRRORA/G230L 

2950 

33752 

20120411 

6955 

cos 

MIRRORA/G230L 

2635 

13651 

20120411 

6955 

cos 

MIRRORA/G230L 

3360 

20101 

20111226 

6848 

WFC3/UVIS a 

F218W 

2228.5 

3000 

20120218 

6902 

WFC3/UVIS 

F275W 

2710.2 

3000 

20120219 

6903 

WFC3/UVIS 

F336W 

3354.8 

3000 

20111224 

6846 

WFC3/UVIS 

F438W 

4326.5 

856 

20111224 

6846 

WFC3/UVIS 

F555W 

5308.1 

856 

20111224 

6846 

WFC3/UVIS 

F625W 

6241.3 

856 

20120219 

6903 

WFC3/UVIS 

F814W 

8029.5 

836 

20120219 

6903 

WFC3/UVIS 

F850LP 

9168.80 

836 

20120219 

6903 

WFC3/UVIS 

F105W 

10552.00 

328 

20120219 

6903 

WFC3/UVIS 

F125W 

12486.00 

328 

20120219 

6903 

WFC3/UVIS 

F160W 

15369.00 

328 


Notes. a For the UVIS images, the UVIS2 aperture was used, so that the SN held was near the edge of chip 2. 
The UVIS exposures were dithered. No post-hash was used for any of this imaging. The IR observations in each 
band were conducted using the SPARS25 sequence with the number of samples equal to 14. 


Table 2 

WFC3 UVIS Photometry 


Target 

F218W 

F275W 

F336W 

F438W 
Vegamag (la 

F555W 

Uncertainty) 

F625W 

F814W 

F850LP 

SN 1993J 

21.657 (0.037) 

21.616(0.017) 

22.250 (0.016) 

22.580 (0.017) 

21.499(0.006)) 

21.373 (0.007)) 

20.885 (0.009)) 

21.393 (0.022) 

A 

23.527 (0.111) 

23.330 (0.041) 

23.414(0.029) 

24.126 (0.039) 

22.681 (0.012) 

21.752 (0.009) 

20.634 (0.007) 

20.321 (0.012) 

B 

23.472 (0.114) 

22.733 (0.029) 

22.657 (0.018) 

23.175 (0.022) 

23.127(0.015) 

23.077 (0.018) 

22.915 (0.026) 

23.043 (0.060) 

C 

23.152(0.097) 

22.354 (0.024) 

22.522 (0.017) 

23.784 (0.030) 

23.767(0.021) 

23.717 (0.026) 

23.713 (0.045) 

23.833 (0.106) 

D 

22.468 (0.057) 

22.071 (0.021) 

22.411 (0.016) 

23.825 (0.031) 

23.868 (0.022) 

23.922 (0.030) 

24.213 (0.062) 

24.506 (0.175) 

E 

23.186 (0.141) 

23.323 (0.041) 

23.238 (0.025) 

24.063 (0.036) 

24.129(0.026) 

24.139 (0.034) 

24.093 (0.057) 

24.143 (0.129) 

F 

24.569 (0.269) 

23.869 (0.058) 

23.840 (0.036) 

25.157 (0.077) 

24.876 (0.045) 

24.897 (0.061) 

24.751 (0.096) 

24.808 (0.220) 

G 

23.495 (0.115) 

23.000 (0.036) 

23.182 (0.025) 

24.372 (0.044) 

24.072 (0.026) 

23.836 (0.029) 

23.357 (0.038) 

23.018(0.062) 

H 

>25.0 

24.070 (0.070) 

24.182 (0.047) 

24.940 (0.065) 

24.890 (0.045) 

24.624 (0.049) 

24.379 (0.078) 

24.007 (0.118) 

I 

>25.0 

24.396 (0.083) 

24.532 (0.057) 

25.954 (0.132) 

25.661 (0.074) 

25.594 (0.098) 

>26.2 

>26.4 

J 

>25.0 

>26.4 

>27.0 

>27.2 

26.569 (0.147) 

25.606(0.102) 

23.825 (0.048) 

23.461 (0.080) 

K 

24.996 (0.346) 

23.394 (0.042) 

23.524 (0.029) 

25.187 (0.109) 

25.033 (0.048) 

24.805 (0.055) 

24.435 (0.071) 

24.959 (0.246) 

L 

>25.0 

>26.4 

>27.0 

25.364 (0.082) 

24.882 (0.042) 

24.522 (0.044) 

24.274 (0.065) 

24.131 (0.128) 

M 

24.615 (0.266) 

23.422 (0.074) 

23.549 (0.029) 

24.740 (0.055) 

24.934 (0.043) 

24.903 (0.058) 

24.561 (0.081) 

25.183 (0.298) 

N 

>25.0 

24.412 (0.080) 

24.360 (0.060) 

25.393 (0.084) 

25.443 (0.062) 

25.455 (0.086) 

24.538 (0.079) 

24.261 (0.142) 


Table 3 

WFC3 IR Photometry 3 


Target 

F105W 

F125W 

F160W 


Vegamag (la Uncertainty) 

SN 1993J/companion 

18.691 (0.005) 

19.530 (0.011) 

19.389 (0.015) 

A 

19.812 (0.005) 

19.276 (0.005) 

18.666 (0.005) 

B 

21.867 (0.018) 

21.192 (0.014) 

20.644 (0.017) 

C 

22.799 (0.038) 

22.385 (0.039) 

21.634 (0.035) 


Notes. a Stars D-N were not detected by Dolphot and have the following 3 a 
upper limits: F105W > 24.8, F125W > 24.2, and F160W > 23.3 mag. 


obtained with the HST Advanced Camera for Surveys’ (ACS) 
High Resolution Camera and downloaded from the Hubble 
Legacy Archive. 11 The comparison of these epochs highlights 
two expanding rings. The expansion of these rings can be 
described similarly to echoes observed in other SNe, remnants, 


11 http://hla.stsci.edu 


and flares (e.g., Rest et al. 2011, 2012 and references therein). 
Owing to a time delay associated with distance, the observer 
witnesses the scattered SN light from different locations on 
large clouds of dust. Here we highlight the detection of these 
echoes, but a complete analysis entails many subtleties and is 
beyond the scope of this paper. 

3. KECK OBSERVATIONS 

On 2013 February 17, we obtained a spectrum of SN 1993J 
with the DEep Imaging Multi-Object Spectrograph (Faber 
et al. 2003) mounted on the 10 m Keck II telescope. It was 
obtained using the 1200/7500 grating, along with a 1" wide slit, 
resulting in a wavelength coverage of 4750-7400 A and a typical 
resolution of 3 A. The slit was aligned along the parallactic angle 
to minimize differential light losses (Filippenko 1982). 

The data were reduced using standard techniques (e.g., Foley 
et al. 2003; Silverman et al. 2012). Routine CCD processing and 
spectrum extraction were completed with IRAF, and the data 
were extracted with the optimal algorithm of Horne (1986). We 
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Figure 2. WFC3 photometry of SN 1993J and the surrounding Stars A-N 
identified in Figure 1 . 

(A color version of this figure is available in the online journal.) 


obtained the wavelength scale from low-order polynomial fits 
to calibration-lamp spectra. Small wavelength shifts were then 
applied to the data after cross-correlating a template sky to an 
extracted night-sky spectrum. Using our own IDL routines, we 
fit a spectrophotometric standard-star spectrum to the data in 
order to flux calibrate the SN and to remove telluric absorption 
lines (Wade & Home 1988; Matheson et al. 2000). 

4. COS OBSERVATIONS 

SN 1993J was observed with the HST Cosmic Origins 
Spectrograph (COS) in 2012 March- April as part of program 
GO- 12531 (PI: A. Filippenko), as summarized in Table 4. For all 
HST COS observations, SN 1993J was centered in the aperture 
using a COS NUV target acquisition (TA) image with MIRROR 
A and the Primary Science Aperture (PSA). Analysis of the TA 
images shows that this strategy was successful and that SN 1993J 
is centered in the aperture during all subsequent observations. 

The COS data consist of far-UV (FUV) and NUV observa- 
tions taken over three visits (11 orbits). The COS NUV obser- 
vations use the G230L grating with central wavelengths 2635, 
2950, and 3360 A. The 2635 and 3360 A central wavelengths 
were observed in a single visit using two FP-POS positions 
(dispersion direction dithers) for each central wavelength to 
correct fixed-pattern noise. The 2950 A central wavelength was 
observed in a second visit using four FP-POS positions. Al- 
though the COS NUV spectroscopic modes disperse light onto 
three stripes on the detector (similar to an echelle grating), the 
third stripe (Stripe C) of the G230L grating is contaminated with 
second-order light and the fluxes from this stripe are unreliable. 
We therefore exclude it from our analysis. The COS FUV data 
were obtained in the final visit using the G140L grating and 
the 1105 A central wavelength. The COS FUV mode typically 



(a) F555Win2002VkA^*AT' 
t'X : .-y* ' ’-Jp feSfciKS 



Figure 3. F555W images highlighting light echoes observed around SN 1993J 
in (a) 2002 and (b) 2012. Both panels are showed in the same orientation as 
Figure 1, with north as up and west to the right. 


Table 4 

HST GO- 12531 COS Spectroscopic Observations 


UT Date 

Epoch 

(days) 

Grating 

Central Wavelength 

(A) 

Exposure 

(S) 

20120331 

6944 

G140L 

1105 

7466 

20120406 

6950 

G230L 

2950 

33752 

20120411 

6955 

G230L 

2635 

13651 

20120411 

6955 

G230L 

3360 

20101 


places the spectrum across two detectors: Segment A and Seg- 
ment B. However, when the 1105 A mode is used. Segment B 
must be turned off to avoid damaging the detector with zeroth- 
order light, so we present only Segment A data. The FUV data 
were taken using three FP-POS positions. 

Figure 1(d) displays the COS field of view. Despite its large 
aperture diameter (2"5), we chose COS for its superior resolution 
and low dark rate compared to Space Telescope Imaging 
Spectrograph (STIS) (see the COS Instrument Handbook). 
Figure 1 highlights that while SN 1993J is the brightest object 
in the aperture, a number of nearby sources are also UV bright. 
Even at this resolution their line-spread functions overlap and 
spatially contaminate the SN spectrum. Spectra were extracted 
using the COS calibration pipeline, CALCOS. 12 The default 
extraction parameters were modified to best isolate the SN 1993J 
flux from other objects in this crowded field. We stress that 
even with HST ’ s superior angular resolution, isolating SN 1993J 
proved difficult owing to both the faint NUV flux and relatively 
low FUV cross-dispersion resolutions (see Figure 4). 

4.1. Defining the Extraction Box 

While the COS low-resolution mode offers superior spectral 
resolution and lower background rates compared to similar 
STIS modes, its large slitless aperture inevitably integrates 
light from sources close to the SN; see Figure 1(d). Figure 4 
plots the cross-dispersion (spatial) profile of all sources within 
the COS aperture, created by summing the two-dimensional 
images in the dispersion direction. The NUV and FUV spectra 
do not fall in a flat line across the detector. For this reason, a 
cross-dispersion profile created from all pixels in the dispersion 
direction will produce broad profiles which do not represent the 
true width of the spectrum on the detector. To balance the spatial 
resolution with the signal-to-noise ratio (S/N), we sum over 
pixels 5000-7000 in the COS/FUV/G140L/1 105 spectrogram 
and pixels 440-520 in the COS/NUV TA image. This figure 
highlights that the default extraction boxes used by the HST 


12 http://www.stsci.edu/hst/cos/pipeline/ 
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Figure 4. Visit 07 TA image zoomed in on the target location (left). The cross-dispersion profiles of the FUV and NUV observations are in green and purple, 
respectively. The extraction box for the spectrum is shown in cyan and the extraction boxes for the background regions are in yellow. In both plots, the locations of the 
stars found in Figure 1 are labeled and marked with horizontal red lines. It is apparent in the FUV cross-dispersion profile that Stars E-I have overlapping line-spread 
functions. 

(A color version of this figure is available in the online journal.) 


Table 5 

Customized Spectral Extraction Parameters 2 


Mode 

93J (F, G, H, I) 
Extraction 
Box Center 
Seg A/Seg B 

93J (F, G, H, I) 
Extraction 
Box Size (Centered 
on Column 1) 

Background 
Extraction 
Location Offsets 
From Spectrum 

Background 
Extraction Box 
Height (Centered 
on Column 3) 

G140L/1105 

496.24/- 

11 

20 

ii 

G230L/2950 

176.04/276.64 

39 

70 

39 

G230L/2635 

177.82/278.19 

39 

70 

39 

G230L/3360 

179.86/280.07 

39 

70 

39 


Note. a All units in pixels. 


pipeline (52 pixels for the NUV and 47 pixel for the FUV) 
integrate point-spread functions (PSFs) from Stars A-I in the 
resulting spectrum. We therefore utilize the NUV TA image to 
optimize our extraction parameters. 

The COS NUV detector has the highest resolution of any 
instrument on HST, with 24 mas pixel 1 and a 3-pixel resolution 
element in both the dispersion and cross-dispersion directions. 
For the NUV spectra, a very narrow (10 pixel) extraction box 
can isolate the SN from the other stars, but this would result in 
significant flux loss from the wings. Since the NUV counts are 
low to begin with ((1.5-3. 5) x 1 0 4 counts s -1 ), the additional 
flux from these wings proves important. Furthermore, the 
FUV/G140L grating has a lower resolution of 80.3 mas pixel' 1 
and a 6 x 10 pixel resolution element (dispersion by cross- 
dispersion). Given the FUV resolution-element size and the 
S/N of the NUV data, we define the extraction box for both 
detectors as ~1 FUV cross-dispersion resolution element of 1 1 
FUV pixels (39 pixels in the NUV; see Figure 4 and Table 5). 
Our resulting spectrum should therefore be considered a blend 


of the SN plus Stars E-I, referred hereafter as “SN+Stars 
E-I.” In Section 5 we will attempt to disentangle the different 
components by combining the photometry and spectroscopy 
with stellar models. 

The dark rate in the NUV detector does not vary with detector 
location. To minimize contamination by other sources, we define 
the background region above Star C and below Star D (see 
Figure 4). Although the lower background region includes Star 
B, placing the window any lower in the NUV results in the 
overlapping of the background region for Stripe A with the 
spectrum for Stripe B. The FUV detector dark rate does vary 
with detector location, but these same background windows 
receive the minimum contamination from other sources. For 
consistency, we therefore define the same background regions 
in the FUV. It should be noted that owing to the faint SN flux, 
the background count rate dominates the total count rate. A 
visual inspection of the background at different extraction box 
heights confirms that the chosen background extraction box 
yields accurate background rates. Table 5 and Figure 4 highlight 
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Figure 5. Ratio of the flux in an 1 1 -pixel-high extraction box with nonstandard 
background regions to the flux in a default extraction as a function of wavelength 
for G140L/1105 sensitivity monitoring observations. The ratio for individual 
observations is in blue, the mean of each observation in red, and a cubic spline 
fit to the mean in green. The mean data values are linearly interpolated in the 
range 1 185-1250 A to avoid fitting the Lya airglow line. 

(A color version of this figure is available in the online journal.) 


the spectrum extraction box center and height, as well as the two 
background region extraction box centers and heights for each 
central wavelength. 

4.2. Extracting with a Nonstandard Box Height 

All extractions are performed using the COS calibration 
software CALCOS 2.20. 1 . When not using the default extraction 
parameters, however, two additional details must be considered. 
First, the flux calibration in CALCOS does not correct for the 
light which falls outside a nonstandard extraction box. Second, 
the default centering of the extraction box assumes that the target 
is centered in the PSA. The location of the PSA on the detector 
is identified with the wavelength calibration line lamp (wavecal) 
spectrum and a known offset between the wavecal aperture 
and the center of the PSA. Since the default extraction box 
is typically large, the centering is only required to be accurate to 
within a few pixels. The SN 1993J extraction requires a much 
smaller extraction box height, so the location of the extraction 
box center must be more precise. 

The SN 1993J COS observations are very faint and blended, 
so they cannot be used to directly calculate either the aperture 
correction or extraction box center. Instead, we utilize the spec- 
trum of bright, well-characterized white dwarfs used to monitor 
the sensitivity of both the COS FUV and NUV detectors over 
time. The bright and well-defined fluxes allow for an accurate 
measurement of the spectrum center and flux corrections for 
small extraction boxes. Furthermore, the large number of regu- 
lar observations throughout COS operations allow us to measure 
the repeatability of these values. 

The CALCOS xld task was run on 14 G140L/1 105 observa- 
tions of WD0947+857 using the extraction parameters detailed 
in Table 5 for SN 1993J. Figure 5 plots the ratio of the FUV flux- 
calibrated spectrum with these custom extraction parameters to 
the CALSPEC model spectrum 1 for each observation. A cubic 
spline is fit to the mean value of the observations as a function 
of wavelength. The uncertainty of the fit at each wavelength 
is calculated by fitting a cubic spline to the standard deviation 
of the observations around the fit. The fit to the flux ratio is 


13 http://www.stsci.edu/hst/observatory/cdbs/calspec.html 
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Figure 6. Plot of spectra extracted for each central wavelength, (a) Individual 
flux spectra, (b) Individual uncertainties, (c) Resulting error bars when flux 
spectra combined using the splice task in the PyRAF STSDAS package, 
(d) Resulting error bars when final spectrum built by choosing only the minimum 
error and associated flux at a given wavelength (i.e., no coadding of spectra). 
The effects of these methods are discussed in Section 5.3. 

(A color version of this figure is available in the online journal.) 
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used to empirically correct the flux of the extracted spectrum 
of SN 1993J. A similar procedure is used for each stripe of the 
NUV detector. 

By default, the xld task extracts the target spectrum from the 
nominal spectrum location combined with the location of the 
aperture on the detector. However, it is possible to ask CALCOS 
to find the spectrum location. This option is used on the white- 
dwarf spectra to determine if a systematic offset exists between 
the nominal spectrum location and the actual spectrum location. 
Offsets of 2 pixels and —3 pixels are found for the G140L/1 105 
and G230L/2635 modes, respectively. 

4.3. Generating the COS Spectrum 

Each observation of SN 1993J is extracted using the param- 
eters listed in Table 5 and the offset location found above. The 
different FP-POS positions of each central wavelength are com- 
bined using the fpavg task in CALCOS to create one spectrum 
for each central wavelength. The flux for each central wave- 
length is then corrected to an infinite aperture and the errors are 
scaled appropriately. Figure 6 plots the flux and observational 
uncertainty for each central wavelength. Spectra from each cen- 
tral wavelength are combined into a single spectrum using the 
splice task in the PyRAF STSDAS package. Prior to combining 
the central wavelengths, regions of low sensitivity and geocoro- 
nal airglow are flagged as poor quality to avoid contaminating 
the final spectrum (see Table 6 for a list of flagged regions). 

Figure 6 shows, however, that the estimates of the obser- 
vational uncertainties associated with the combined spectrum 
are not well-behaved at wavelengths where the FUV and NUV 
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Table 6 

Wavelength Ranges Excluded from Final COS Spectrum 


Mode 

Segment/ Stripe 

Flagged 

(A) 

Reason 

G140L/1 105 

FUVA 

1198-1202.5 

Geocoronal N i 

G140L/1105 

FUVA 

1210.5-1220.5 

Geocoronal Lya 

G140L/1 105 

FUVA 

1300.5-1307.4 

Geocoronal O i 

G140L/1 105 

FUVA 

1354.3-1359 

Geocoronal O i 

G140L/1105 

FUVA 

A > 1665 

Low sensitivity 

G230L/2635 

NUVA 

X < 1665 

Low sensitivity 

G230L/2950 

NUVB 

X > 3000 

Low sensitivity, second-order light 

G230L/3360 

NUVB 

X > 3000 

Low sensitivity, second-order light 



Figure 7. COS FUV and NUV observations combined into a single spectrum 
with regions of poor data quality excluded. The spectrum is in black with 
uncertainties in gray. The spectrum has been boxcar smoothed by 70 pixels. 


spectra overlap. A detailed analysis in Section 5.3, below, fur- 
ther shows that over large regions of the spectrum, observational 
uncertainties are overestimated by varying factors depending on 
which central-wavelength grating is used. Instead, we find that 
better-behaved observational uncertainties can be obtained by 
not coadding spectra at wavelengths where the FUV and NUV 
spectra overlap. We choose to extract only the minimum obser- 
vational uncertainty and associated flux at these wavelengths. 
The excluded regions are listed in Table 6. Figure 6 also plots 
the resulting uncertainties for this technique. In Section 5.3, we 
further show that these results are substantially more uniform 
over the entire spectral range, and the overall distribution is very 
close to normally distributed. 

Figure 7 plots the final COS spectrum, which should be 
considered a sum of flux from the SN+Stars E-I. These are the 
same sources that must have been included in the Keck/LRIS 
spectrum presented by Maund et al. (2004). Unlabeled fainter 
stars are apparent in the WFC3/F336W image in Figure 1, but 
are not detected in the COS NUV image. We thus conclude 
that these stars do not contribute significantly to the extracted 
COS spectrum. 

4.4. Spectral Line Identifications and Profiles 

Figure 8 shows a compilation of the most prominent UV lines 
in the COS spectrum. We also include several optical lines from 
the 2013 February 17 Keck spectrum (see Section 3). In this 
latter spectrum the Ha and [Om] a/,4959, 5007 emission lines 
are strongest. Even the blended [O in] doublet shows consider- 
able structure, which reflects real fluctuations in the emissivity 
as a function of the velocity. The reality of these fluctuations 



Velocity (1 0 3 km s ') 


Figure 8. Compilation of UV and optical emission-line profiles. The flux of each 
line is normalized within the velocity range ±10,000 km s -1 . The positions of 
the Sin 7.1526 and Civ interstellar absorptions are marked in the Civ panel. 
Dashed lines indicate [Om] line peaks at 3800 km s -1 and +3700 km s -1 . 
Dashed lines also mark the constraints on the ejecta velocity given by the Ha 
line at 6000 km s -1 and ~+9000 km s -1 . 

(A color version of this figure is available in the online journal.) 

is demonstrated by the same two peaks at —3472 km s -1 and 
—4272 km s' 1 in both the 14959 and 15007 components of 
the [O in] doublet. The [O hi] line shows an additional asym- 
metry favoring the blue side, which agrees with observations of 
SN 1993J by Milisavljevic et al. (2012). 

Several other lines show structures similar to that of [Om]. 
In particular, the peak at 3800 km s 1 is clearly present in 
the Ha, Lya, [Oi], Cm], and Civ lines. Note the interstellar 
absorption of the Sin 11526 line, which distorts the Civ peak. 
The N n], N hi], and Mg ii lines have a lower S/N, and for these 
lines the ~— 3800 km s _I peak is the most significant extension 
above the noise. Dashed lines in Figure 8 highlight the strongest 
peaks at both ~— 3800 km s -1 and ~+3700 km s” 1 . These 
velocities also correspond to the flat part of the box-like line 
profiles previously identified by Fransson et al. (2005). 

The asymmetry is less pronounced in the Ha line compared to 
the [O m] line, although the fluctuations are present. These lines 
likely arise from different components, with the Ha line coming 
mainly from the dense, cool shell behind the reverse shock, while 
the [O hi] and other high-ionization lines come from the ionized 
outer parts of the unshocked ejecta (e.g., Chevalier & Fransson 
2003). The full extension of the lines, marking the maximum 
ejecta velocity, is somewhat difficult to determine because of the 
gradual transition to the continuum. The red wing of the Ha line 
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Figure 9. UV spectrum of SN 1993J from days 2579-2585 and 6944-6955. 
The most prominent emission lines have been marked. The top spectrum has 
been shifted up by 1 x 10 16 erg cm” 2 s” 1 A” 1 for clarity. 

(A color version of this figure is available in the online journal.) 


Table 7 

Wavelength Ranges Excluded from Final STIS Spectrum 


Mode 

Flagged 

Reason 


(A) 


G140L 

7 < 1141 

Low sensitivity 

G140L 

1206-1220 

Geocoronal Lya 

G230L 

7 < 1700 

Low sensitivity 


extends to ~+9000 km s _1 , although the maximum velocity may 
be influenced by the [N ii] 7.6583 line. The blue wing is blended 
with the [O i] 7,7.6300, 6364 lines, which allows us to constrain 
an ejecta velocity of at least ~6500 km s 1 on this side. Figure 8 
also highlights these velocities with dashed lines. Compared to 
high-ionization lines (i.e., the Civ and [Oni]), Ha has a larger 
width on both the red and blue sides, with a maximum velocity 
of 5000-6000 km s _1 . As discussed above, these velocities may 
arise from different regions where these lines are formed. In 
particular, the high-ionization lines are likely to be formed by 
photoionization in the processed gas in the unshocked core, and 
are consequently expected to have lower velocity that the Ha 
line. 

5. ANALYSIS 

Maund et al. (2004) detected a number of absorption lines in 
the NUV portion of their Keck/LRIS spectrum of SN 1993J. 
Models indicate that these lines are most consistent with those of 
a B2 la star at the position of SN 1993J. Maund & Smartt (2009) 
later found a NUV excess in an HST SED consistent with the 
combined flux from a B2 la star and SN 1993J, although the SN 
still dominated the SED. Maund & Smartt (2009) predicted that 
by 2012, the SN would fade sufficiently to directly measure 
the properties of the putative remaining hot B companion 
star. Here we disentangle the contributing sources in the 2012 
COS spectroscopy. 

5.1. Defining a Supernova Template 

Since the SN and putative companion are spatially coincident, 
we must first define a template representative of the SN 



Figure 10. COS UV spectrum fit with only the year-2000 STIS spectrum scaled 
by a factor of 0.3 1 . This scaled STIS fails most significantly in fitting the rising 
FUV continuum. Subsequent figures describe additional components. 

(A color version of this figure is available in the online journal.) 


contribution to Ihe overall spectrum. Prior to 2012, the most 
recent UV observations of SN 1993J were taken using the STIS 
in 2000 April with the G140L and G230L gratings (GO-8243; 
PI: R. Kirshner; Fransson et al. 2005). This spectrum serves 
as a useful reference template; we can assume that by this 
epoch the SN was quite evolved, and we expect little change 
in the spectrum besides continued fading. The validity of this 
assumption is highlighted by the slow UV spectral evolution 
in 1995-2000 shown by Fransson et al. (2005, their Figure 1) 
and the slow optical spectral evolution in 2000-2009 shown 
by Milisavljevic et al. (2012, their Figures 10 and 11). While 
some line strengths do vary, the underlying continuum shape 
and line profiles of the strongest emission lines do not change 
significantly, which is most important for this study. 

For the 2000 STIS observations, the SN flux dominated 
over all neighboring stars, so contamination is not an issue. 
We therefore use the default pipeline product, just as Fransson 
et al. (2005) did. These observations are again combined using 
the splice task. Like the COS spectrum, wavelengths of low 
sensitivity and geocoronal airglow are flagged prior to the 
combination. A list of excluded wavelength ranges can be found 
in Table 7. 

Figure 9 compares the 2000 STIS and 2012 COS spectra. 
Although the S/N is considerably lower in the 2012 spectrum, 
most of the lines identified in 2000 are also observed in 2012. As 
in 2000, the strongest line is Lya, but it is severely affected by 
the strong geocoronal emission line, the geocoronal subtraction, 
and the damping wings of the interstellar Lya absorption. We 
therefore do not include this line in our analysis, although we do 
note that it appears asymmetric with a strong blue wing. Both 
the Civ 7.7.1548, 1551 and Chi] 71909 lines are still strong in 
the 2012 spectrum. The Nil] 7.7.2140,2144 doublet is weaker, 
but still above the noise. The Nm] 7.7.1747-1754 multiplet is 
in a noisy region of the spectrum, but the blue peak is above 
the noise. In Fransson et al. (2005), the feature at 2335 A was 
identified with Cii] 77.2323, 2328, Oiii] 772321, 2331, and 
Sin] 72335, and it is also well defined in the 2012 spectrum. 
The peaks at ~2460 A and ~2640 A are probably produced by 
Fe ii resonance lines. The most dramatic change is the strength 
of the Mg ii 772796, 2803 doublet, which was the strongest UV 
line in 2000 but is only barely above the noise in 2012. 
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Figure 11. COS UV spectrum is better fit by including stellar models for Stars E, F, G, and H along with the year-2000 STIS spectrum scaled by a factor of 0.28. 
Nonetheless, the FUV continuum excess remains unaccounted for. Top: the entire / 1ST data set, including COS spectra and WFC3 photometry. The longer-wavelength 
photometry constrains the stellar models used for Stars E-H. The excess flux for SN 1993J at these longer wavelengths is attributed almost entirely to the SN component 
and not the putative companion; we therefore do not fit these data. Bottom: a closer look at the UV spectrum and best-fit model. 

(A color version of this figure is available in the online journal. 1 


5.2. Modeling the COS Spectrum 

Figure 7 plots the combined COS NUV and FUV spectrum. 
We again stress that the spectrum should be considered a sum 
of flux from the SN+Stars E-I (see Section 4). Since the optical 
wavelengths are dominated by the fading SN itself, we aim here 
to fit the combined NUV and FUV spectrum (i.e., k < 3000 A). 
This wavelength range is not sensitive to the same stellar 
absorption lines observed by Maund et al. (2004) at >3600 A 
(see spectral models below). To model the data, we consider 
all possible components with the intent of minimizing x 2 (see 
Section 5.3). 

Fading SN. As described in Section 5.1, we adopt the UV 
spectrum obtained with HST /STIS in 2000 (Fransson et al. 
2005) as the template of the fading SN 1993J. We use the IDL 
MPFIT nonlinear least-squares curve fitting function (Markwardt 
2009) to find the best fit by varying a single multiplicative 
scale factor across the entire spectrum. Figure 10 shows that a 


best fit can be achieved with a scale factor of 0.31 times the 
original STIS flux in 2000. At longer wavelengths (the NUV), 
the continuum levels match well. The lines, however, vary in 
their strength and do not appear to be well represented by a 
single scale factor. At shorter wavelengths, an FUV continuum 
excess exists which can be best fit by additional components 
described below. 

Stars E-H. Figure 11 plots the WFC3 photometry for 
SN 1993J and Stars E-H (Star I was too faint to be detected 
by Dolphot). The COS spectrum blends Stars E-I together 
with the SN, but the WFC3 photometry isolates the flux from 
the individual stars. For this reason, the SN 1993J photometry 
does not match the spectrum. The combined flux from all the 
stars, however, is consistent with the spectrum to within the 
uncertainties. 

We fit Stars E-H using Castelli and Kurucz models (Castelli 
& Kurucz 2004) convolved with the WFC3 throughput for each 
filter using IRAF’s synphot.calcspec package. For simplicity, 
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Figure 12. Same as Figure 11, except with an additional component representing a B2 la T = 20,000 K stellar model normalized to V = 22 mag (orange). While the 
overall fit is better than in Figure 1 1, the B2 la model is brighter than the WFC3 photometry for SN 1993J. 

(A color version of this figure is available in the online journal.) 


we assume the same extinction as the SN (see below). The model 
chosen for each star is the model for which x 2 is minimized. 
These stellar models provide reliable estimates of the FUV 
fluxes despite a lack of photometry at these wavelengths. 
Figure 11 shows the resulting best fits. We then rescale the 
STIS spectrum to account for the remaining flux. Combined 
with the flux from Stars E-H, an STIS scale factor of ~0.28 
results in a best fit, but there is still an FUV excess implying the 
existence of another component. 

Additional B2 star. The contribution of FUV flux from Stars 
E-H is minimal, and even with the addition of these stars, an 
FUV excess is still apparent. We therefore consider the presence 
of a binary companion, previously suggested to be a B2 la star 
(Maund et al. 2004; Maund & Smartt 2009). We generate a 
B2 la spectrum using the Castelli and Kurucz model in IRAF’s 
synphot.calcspec package. We use a 20,000 K star with a surface 
gravity log(g/cm s -2 ) = 3.5 and metallicity [M/H] = 0.0. 
Following Maund et al. (2004), we renormalize to the observed 
V-band magnitude ( VVegamag = 22.0 accounting for an extinction 
of Ay = 0.52 mag to SN 1993J). Again, we use the IDL MPFIT 
function to find the best fit by varying only the multiplicative 
scale factor for the STIS spectrum. 


Figure 12 shows an improved fit with an STIS scale factor 
of ~0.12. This model, however, is not ideal. The combined 
flux from the scaled STIS spectrum and B-star model (orange) 
should be consistent with the observed WFC3 photometry, 
but in this case the B-star model is too bright. To optimize 
our fit, we consider other stellar types and luminosity classes 
ranging from 10,000 K to 31,000 K, as well as normalization 
magnitudes (although we do not vary the metallicity or surface 
gravity values). We minimize / 2 (e.g., maximize the probability 
function) with a temperature of 24,000 K renormalized to 
Vvegamag = 22.9 mag (see Figure 13). Given the noise associated 
with the spectrum, the data are somewhat degenerate with 
several models that include B stars in the range 19,000-27,000 K 
(B3-B1). We thus infer the presence of a spectral component 
consistent with continuum from a hot B star. 

5.3. Minimizing Reduced x 2 

To effectively utilize the reduced x 2 requires a detailed 
understanding of the observational uncertainties. CALCOS, 
the COS calibration pipeline, generates both a flux-calibrated 
one-dimensional spectrum and an associated uncertainty. This 
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Figure 13. Same as Figure 12, except with a hotter B1 T = 24,000 K stellar model 
model while remaining consistent with the WFC3 photometry for SN 1993J. 

(A color version of this figure is available in the online journal.) 

uncertainty is calculated using the formula given by Gehrels 
(1986; see COS Data Handbook; Massa & York 2013). Gehrels’ 
formula, however, offers only an upper limit to the uncertainty 
regardless of the number of source counts (for large counts the 
formula approaches the typical *Jn error approximation). 

Figure 6 shows two methods used for generating our spectrum 
of SN 1993J. In the first, we combine all spectra from each cen- 
tral wavelength into a single spectrum. In the second, we do not 
coadd spectra at wavelengths where the FUV and NUV spectra 
overlap. We choose to extract only the minimum observational 
uncertainty and associated flux at these wavelengths. For each 
case, the residual errors can be characterized by the distribution 
of x = (Observed(i) — M o d c I (/) ) / U n c e r t a i n t y ( / ) , for which a 
normal distribution should be characterized by a Gaussian with 
a standard deviation equal to 1 . For the two techniques described 
here. Figures 14(a) and (b) show the values / as a function of 
wavelength and the resulting histogram. 

Figure 14(a) varies significantly, while Figure 14(b) shows a 
more normal distribution. Still, the uncertainties appear too large 
in the FUV (i.e., / < 1). To compensate, we decrease the size of 
all uncertainties in spectral elements with a wavelength shorter 


normalized to V = 22.9 mag (orange). The value of / 2 was minimized for this 


than 1660 A by a factor of 2.9, and all spectral elements with a 
wavelength longer than 1660 A by a factor of 1.4. Figure 14(c) 
shows the results, which is nearly a normal distribution. We 
adopt these uncertainties throughout the modeling in this paper. 

Despite the improved accuracy of the uncertainties, the 
resulting reduced / 2 values cannot be considered meaningful 
on an absolute scale because the degrees of freedom are not 
well defined. We therefore calculate the probability (Q) that the 
resulting / 2 value is due to chance and normalize relative to the 
sum of the resulting probability distribution (see Figure 15). We 
consider models summing to 68% of the combined probability 
(e.g., 19,000-27,000 K) to be the likely range of models that 
can fit the COS spectrum. 

5.4. Other Sources of FUV Flux 

While we do find, like Maund et al. (2004), that the surround- 
ing stars contribute little to the overall flux, we also consider the 
possibility that this hot B-star emission may not arise from the 
supposed binary companion. First, we cannot definitively rule 
out line-of-sight coincidences. Moreover, the high-resolution 
TA image contours in Figure 4 reveal the presence of an 
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Figure 14. For each method used to generate the final COS spectrum, 
the etrors are represented by the distribution of '/ = (Observed(Z) — 
Model(i))/Uncertainty(i), for which a normal distribution should be charac- 
terized by a Gaussian with a standard deviation equal to one. (left) / vs. wave- 
length. (right) Histogram of the x distribution, overplotted with a Gaussian with 
a standard deviation equal to one. (a) Results from using the splice task in the 
PyRAF STSDAS package, (b) Results when selecting only the minimum error 
and associated flux at a given wavelength (i.e., no coadding), (c) Same as (b), 
but the errors are manually corrected to achieve a normal distribution with a 
standard deviation close to one. Ultimately, (c) illustrates the method used to 
generate the COS spectrum used for our modeling in Section 5.2. 

(A color version of this figure is available in the online journal.) 



Figure 15. Probability ( Q ) that the best-fitting model / 2 value is due to chance, 
normalized relative to the sum of the resulting probability distribution. Models 
summing to 68% of the combined probability range from 19,000 K to 27,000 K. 


unresolved UV source just south of the SN position. We also see 
a hint of this feature in the WFC/F336W filter in Figure 1 . At a 
distance of 3.6Mpc, this angular separation of >0'.' 1 corresponds 
to >1.5 pc, thereby ruling out the possibility of this source 


being the putative companion of the progenitor. Furthermore, the 
flux of this unresolved source is low relative to the SN+putative 
companion. The presence of this emission, however, suggests 
the possibility that the local progenitor environment is crowded. 
If this is true, the B2-star detection within the HST PSF does 
not necessarily have to be the binary companion. Future HST 
imaging may be able to probe the progenitor environment in 
more detail once the SN sufficiently fades. 

6. CONCLUSION AND FUTURE WORK 

We present COS spectra and WFC3 photometry of SN 1993J 
nearly 20 years post-explosion. The SN emission has faded 
sufficiently that we can test predictions of the presence of the 
progenitor’s companion. The rise in the NUV, previously ex- 
plained by the presence of the hot star, is still consistent with 
the scaled spectrum of SN 1993J obtained by HST/STIS in the 
year 2000. The continued rise in the FU V, however, cannot be ex- 
plained by this model, even with flux contributions from neigh- 
boring stars that contaminate the spectrum. The addition of a hot 
B-star (B3-B1) spectrum most accurately accounts for this FUV 
excess. 

The properties of this companion can be used to directly 
address several questions. The temperature, surface gravity, and 
bolometric luminosity can be used to constrain the original 
stellar mass (e.g., Claeys et al. 2011). The B2 star binary 
companion scenario is also interesting when considering the 
SN 1993J circumstellar medium (CSM), which has a flattened 
or disk-like geometry (Matheson et al. 2000), and derived 
progenitor mass-loss rates of 10 4 — 10 — 5 M Q yr " 1 for a wind 
velocity of 10 km s _1 (e.g., Van Dyk et al. 1994; Fransson 
et al. 1996; Mioduszewski et al. 2001; Immler et al. 2001; 
Fransson et al. 2005; Chandra et al. 2009). B2 la stars typically 
have mass-loss rates of (0.5-1) x 10 -6 M Q yr -1 and a wind 
velocity of ~500 km s -1 (Crowther et al. 2006). The fact that 
the CSM density of SN 1993J is three orders of magnitudes 
higher than that of typical B2 stars implies that the SN 1993J 
CSM was formed by the primary star that exploded, or during 
Roche-lobe overflow (RLOF) prior to the explosion. Such a 
system is relevant in the context of several nearby B supergiants 
surrounded by ring-like nebulae, including the B3 1 progenitor 
to SN 1987A (e.g., Burrows et al. 1995; Smith et al. 2007), Sher 
25 (a B1.5 la star; Brandner et al. 1997; Smith et al. 2007), 
HD 168625 (Smith 2007), and SBW1 (a B1.5 la star; Smith 
et al. 2007, 2013). Many of these systems have previously been 
considered products of RLOF from a companion, products of 
mergers, mass loss from a rapidly rotating star, or a stripped- 
envelope mass loser that has already exploded as an SN Ibc or 
lib, such as SN 1993J. None of these possibilities has yet been 
confirmed. 

While this paper disentangles the NUV contributions of 
nearby stars and the underlying SN, the FUV remains relatively 
unconstrained. Now that the SN flux has faded sufficiently, 
future FUV /NUV imaging of SN 1993J is necessary to confirm 
the presence of a hot B star and constrain the FUV contribution 
from these sources. 
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ESA Hubble Space Telescope, obtained from the Space Tele- 
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ciation of Universities for Research in Astronomy (AURA), 
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